function neutral_stab(k)
addpath ../

p = params;

%-----
% neutral stability curves

make_path(2);
% hold on;
% plot([0 5], [p.Ma, p.Ma], 'k-.');
ylim([1e4 1e6]);

